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ABSTRACT 

We investigate the radiative transport of dust in primordial galaxies in the presence of the UV radiation 
field from the first metal-free stars. We find that dust created in the first supernova (SN) explosions can 
be driven through the interior of the SN remnant to accumulate in the SN shells, where second-generation 
stars may form from compressed cooling gas. This scenario requires metal-free stars to form continuously 
over timescales of up to 10 Myr, consistent with recent estimates. Silicate and graphite grains, as well as 
iron-bearing magnetites, are transported to the shells for reasonable parameter assumptions, but their 
relative yields from primordial SNe is an important factor in the resulting abundance ratios. We compare 
the results of segregated grain transport with the current nucleosynthetic data on extremely metal-poor 
Galactic halo stars. Fossil signatures of this process may already have been detected in those iron-poor 
stars with enhanced carbon and silicate elements such as magnesium, silicon and oxygen. We discuss the 
implications of our results for the transition from first- to second-generation star formation in primordial 
galaxies, and the role played by the radiative transport of dust in this process. 

Subject headings: cosmology: theory — dust — galaxies: high-redshift — stars:abundances — stars: 
Population II supernova remnant 



1. INTRODUCTION 

The radiative transport of dust and grain history in 
galaxies have been studied by many authors, generally in 
the context of our Galaxy (e.g., Draine & Salpeter 1979a; 
Dwek & Scalo 1980; Seab & Shull 1983). Although early 
calculations (Pecker 1972) of this problem indicated that 
dust transport was difficult to accomplish, Ferrara et al. 
(1991) found that it was possible to transport dust grains 
by radiation pressure to large distances from the plane of 
the disk in a relatively short time (~ 10 8 yr). In addi- 
tion, dust destruction mechanisms are important to con- 
sider (Shustov & Vibe 1995). 

In general, smaller dust grains are destroyed at a faster 
rate while they are being transported, and they typi- 
cally do not survive. Grain destruction in the interstel- 
lar medium (ISM) occurs through thermal sputtering in 
hot gas as well as non-thermal processes (sputtering, shat- 
tering, vaporization in grain-grain collisions). The non- 
thermal processes occur primarily when the grains are 
swept up by shock waves (Shull 1978; Jones et al. 1994). 
In these situations, the original power-law grain size dis- 
tribution will be modified by various processes of grain 
destruction (Seab & Shull 1983; Jones et al. 1996). Grains 
can also be accelerated through the gas by radiation fields. 
However, at a given radiation pressure, it is more difficult 
to drive large grains owing to their lower ratio of area to 
mass. The final outcome of the radiative transport of dust 
therefore depends on the size and composition of the grains 
as well as the details of the radiation field. 

More recently, numerical simulations in a cosmological 



context have studied the radiative ejection of dust grains 
to the intergalactic medium (IGM) from primordial halos 
in the wake of the first generations of supernovae (SNe) 
(Aguirre et al. 2001; Bianchi & Ferrara 2005). Testing the 
predictions of such theories of cosmological dust transport 
has usually been in a macroscopic context, using upper 
limits from the cosmic microwave background to constrain 
IGM dust distortion (Loeb & Haiman 1997; Ferrara et al. 
1999) or to offer alternative explanations to a cosmological 
constant scenario of the data on Type la SNe at redshifts, 
z < 2 (Aguirre 1999). 

In this paper, we focus on more local effects of dust 
transport within the host galaxy related to second- 
generation stars that form coevally with or subsequent to 
primordial stars. We examine, under the conditions in 
typical primordial halos, whether dust can be driven effec- 
tively by the radiative pressure of metal-free stars. If so, 
can dust grains accumulate in the cooling shells of first- 
generation SN remnants (SNRs), becoming incorporated 
eventually into new star-forming sites? Recent work by 
Mackey et al. (2003) and Salvaterra et al. (2004) suggest 
that such second-generation star formation is highly likely, 
forming from metal-poor or even metal-free fragmenting 
gas in SN shells in early galaxies. The latter authors show 
that instabilities in the SN shell that could lead to gas 
fragmentation and to the eventual formation of low-mass 
second-generation stars set in at about 1-50 Myr after the 
SN explosion. This range holds for the condition that the 
SN shells do not sweep the baryons out of a galaxy of 
virial temperature 10 4 K. A timescale requirement of ~ 
a few Myr to 10 s yr is entirely consistent with the dura- 
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tion of metal-free star formation calculated by both semi- 
analytic (Tumlinson et al. 2004) and numerical methods 
(Wada & Venkatesan 2003; Bromm et al. 2003). These 
involve estimates of the timescales over which metals from 
the first SNe can enrich the gas in either the host galaxy or 
neighboring galaxies. We motivate this hypothesis by con- 
sidering the currently available element abundance ratios 
of extremely metal-poor (EMP) stars in the Galactic halo. 
The elements comprising the dominant dust compounds 
from the first SNe, including C, Si, O, and Fe, in addition 
provide the most effective cooling channels for primordial 
star-forming gas (Bromm & Loeb 2003; Santoro & Shull 
2005). 

The paper is organized as follows. In § 2, we present the 
current nucleosynthetic data on EMP stars to motivate the 
consideration of a dust-transport scenario. In § 3, we de- 
scribe the formalism and assumptions of the model used 
to solve for dust grain transport, the results of which are 
presented in § 4. We discuss implications of our findings 
and conclude in § 5. 

2. NUCLEOSYNTHETIC DATA ON EMP STARS 

We begin by using the current data on EMP stars rel- 
evant for the problem in this paper. We highlight the 
trends of carbon, silicon, oxygen, nitrogen, magnesium, 
and aluminium as a function of iron abundance below 
[Fe/H] ~ —2.8. These elements, with the exception of 
nitrogen, are the most relevant for dust creation from the 
metals from the SNe of metal-free stars (Schneider et al. 
2004; Todini & Ferrara 2001). We define [Fe/H] as the ra- 
tio of the measured column densities of Fe/H to the solar 
ratio, (Fe/H) Q = 4.68 x 10~ 5 . 

The nature of the primordial stellar initial mass func- 
tion (IMF) is currently of great interest and debate. Some 
recent theoretical studies indicate that this IMF may have 
been top-heavy (Abel et al. 2000; Bromm et al. 2002), 
leading predominantly to stellar masses > 100 Mq up to a 
critical gas metallicity of Z CT w 10~ 4±1 Zq (Bromm et al. 
2001a; Schneider et al. 2002), above which a present-day 
IMF occurs. However, other detailed studies of the cur- 
rent data on reionization, high-z star formation and the 
metal abundance ratios in the IGM and EMP stars suggest 
that the primordial IMF, rather than being biased towards 
high masses, may merely lack low-mass stars (Venkatesan 
& Truran 2003; Tumlinson et al. 2004; Daigne et al. 2004; 
Qian & Wasserburg 2005). It is certainly possible that 
both IMFs were coeval in the past, given the right combi- 
nation of conditions, as pointed out by these papers. In ad- 
dition, the transition metallicity Z cr may be significantly 
higher at low densities and could vary with metal species 
(Bromm & Loeb 2003; Santoro & Shull 2005). These two 
IMFs, > 100 M Q and ~ 10-100 M , represent two pos- 
sibilities in the definition of a top-heavy IMF, where the 
IMF's lower or upper mass limit is increased. A third op- 
tion would be to flatten the slope or alter the shape of the 
stellar IMF. 

Regardless of such IMF issues, we expect the first gener- 
ations of stars to form from metal-free gas. Their compo- 
sition heavily influences their structure and properties, as 
they rely predominantly on the p-p chain initially than on 
the more efficient CNO cycle for their thermonuclear fuel 
source (Tumlinson & Shull 2000). Consequently, metal- 



free stars are hotter and emit significantly harder ionizing 
radiation relative to their finite- Z counterparts (Bromm 
et al. 2001b; Tumlinson et al. 2003; Schaerer 2002). This 
will play an important role for the dust transport problem 
here. 

Stars of masses ~ 10-100 Mq end their lives as the more 
familiar Type II SNe, leaving behind neutron stars and 
black holes, whereas metal- free stars in the mass range 
~ 140 — 260M© are thought to disrupt themselves entirely 
as pair-instability SNe (PISNe). We do not concern our- 
selves here with the exact type of SN in primordial halos, 
and require only that the parent metal-free stellar cluster 
has effective luminosities of ~ 1O 6 L (Tumlinson et al. 
2003; Bromm et al. 2001b; Schaerer 2002), which is easily 
achieved by a few tens of stars in a Salpeter IMF in the 
1-100 Mq range or a single star of mass > 100 Mq. For 
the purposes of this paper, we require the former case, i.e., 
a stellar cluster, so that there is a hard photon source for 
at least a few to 10 Myr after the initial SN explosion. We 
also require that the SN kinetic energy is at least 10 51 erg, 
which is true for "normal" Type II SNe. We discuss the 
case of hypernovae (HNe) below; the SN explosion energies 
for HNe and PISNe are thought to lie in the range 10 51 - 
10 53 erg (Umeda & Nomoto 2003; Heger & Woosley 2002). 
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Fig. 1. — The observed abundances of C, N, O (upper panel) and 
Mg, Al and Si (lower panel) relative to Fe as a function of [Fe/H] in 
EMP halo stars. See text for discussion. 

In Figure 1, we show the measured EMP stellar abun- 
dances of C, N, O, Mg, Al, and Si relative to Fe (normal- 
ized to the solar ratio), as a function of [Fe/H]. We limit 
the data to those EMP stars with [Fe/H] < -2.8, to be 
consistent with the upper limit to the transition metallic- 
ity factoring in typical errors in the data of order 0.1-0.2 
dex. This is also close to the upper limit in metallicity de- 
rived by Salvaterra et al. (2004) for second-generation stars 
forming in the SN-induced scenario. We take the data of 
EMP giants from the study by Cayrel et al. (2004) , and of 
EMP dwarfs and other stars from Cohen et al. (2004), Aoki 
et al. (2004), and Tumlinson et al. (2004) and references 
therein. For the two most iron-poor stars, HE-0107-5240 
and the newly discovered HE-1327-2326, we use the pub- 
lished abundances by Christlieb et al. (2004) and Frebel 
et al. (2005). 

We note some relevant points of interest in Figure 1. 
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First, although the two most iron-poor stars have greatly 
enhanced values of C, N and O, there is no overall clear 
trend of [C/Fe] with [Fe/H]. We also observe that those 
EMP stars with highly enhanced C (so-called C-enhanced 
EMPs) 4 also have strong N overabundances relative to so- 
lar, sometimes even exceeding that of C. Oxygen, on the 
other hand, is usually enhanced, sometimes up to 10 2 -10 4 
times solar at the lowest Fe-metallicities. We note also 
the consistent enhancements of Mg and Si of a few to 100 
times solar, and the relatively flat mild underabundance 
of [Al/Fc] with [Fe/H]. 

Several models have been proposed to explain the diver- 
sity of element abundances in low-mass EMP stars. Some 
interpret the abundances as reflecting the star's gas for- 
mation conditions, while others invoke processes subse- 
quent to its creation. The presence of Mg and heavier ele- 
ments would argue for the former class of models in prin- 
ciple. Amongst these, the HN model (Umeda & Nomoto 
2003) is currently the most successful at reproducing the 
EMP data, requiring a new class of SNe at primordial 
metallicities that have enhanced SN explosion energies of 
10 51 — 10 53 erg and enhanced CNO element generation rel- 
ative to Fe at low [Fe/H]. HN models can reproduce most 
EMP element abundances well except N and Na, but they 
require significant fine-tuning for each EMP star, with the 
appropriate ratio of mixing, fallback, asymmetry and/or 
jets in the SN model. The post-formation scenarios in- 
clude the dredge-up of elements during the star's giant 
branch phase, accretion from the interstellar medium or a 
companion star, and the rotation of the parent star - all 
of these could enhance C, N, O and sometimes Mg given 
enough time. The joint trend of enhanced Mg, Si and O 
at low [Fe/H] however remains unexplained. The strong 
variations between individual EMPs amongst the elements 
shown here and those related to the Fe-peak and r- and 
s-process elements imply that any successful scenario may 
require considerable tailoring for each star. 

Although N is not an clement of much relevance to dust, 
we include it here, as it is usually a good tracer of C and 
O, owing to its origin in the CNO cycle, and hence of 
the AGB scenario. We note that we have excluded stars 
known or thought to be in binaries, so we do not neces- 
sarily expect to see the signature from AGBs in binaries 
in Figure 1. Unfortunately, the data on N in this range of 
[Fe/H] is fairly scant, and no firm trend can be stated. As 
noted by Plez & Cohen (2005), EMP stars with enhanced 
C usually have enhanced N as well, but the reverse is not 
always true. This is seen in Figure 1. It is possible that 
N depletion can occur through the formation of ammonia 
(NH3) ices as mantles on grains in dense clouds, but these 
likely are sputtered away immediately. 

In contrast, although not shown here, it is worth noting 
that for the few stars with measured CNO that have [Fe/H] 
> —3, [C/Fe] is mildly underabundant, whereas [N/Fe] is 
mildly overabundant, an anticorrclation noted by Briley 
et al. (2004) in a study of Galactic globular clusters. The 
varied behavior of N/C and the lack of a trend with [Fc/H] 
in EMP stars is noted in Plez & Cohen (2005), who point 
out that the SN origin is attractive owing to the additional 



enhancement of Na, Mg, and especially O in these stars. 
The increasing enhancement of N accompanied by a de- 
cline of C as the star's evolutionary state advances would 
argue for self-enrichment from dredge-ups and other pro- 
cesses in the giant branch phase, but this is not observed. 

In summary, the data on EMP stars indicate that Si, 
Mg, and O, and to a lesser degree C, are generally en- 
hanced in the iron-poor stars. These elements form the 
composition of the dominant dust compounds created in 
SNe from metal-free stars of masses ~ 10-260 Mq (Todini 
& Ferrara 2001; Nozawa et al. 2003; Schneider et al. 2004), 
which include graphites, silicates (enstatites (MgSiOs), 
forsterites (Mg2SiOzt), and Si02), and iron-bearing mag- 
netites. We therefore proceed to investigate whether these 
elements can be selectively transported in a hot radiation 
field, decoupled from the background SN metals in the gas 
phase. 

3. CREATION AND TRANSPORT OF DUST GRAINS 

In this section, we describe the scenario in which dust is 
created in primordial SNe, followed by a summary of the 
assumptions and equations we use to solve for the trans- 
port of dust grains. 

3.1. Dust Formation 

As discussed earlier, the duration of metal-free star for- 
mation is thought to be ~ 10-100 Myr. The lower limit 
of a few-10 Myr has been derived using semianalytic cal- 
culations of halo self-enrichment from SN ejecta (Tumlin- 
son et al. 2004), and from 3D gas hydrodynamic simula- 
tions of the chemodynamical feedback from the first SNe 
on parsec scales in the ISM of primordial galaxies (Wada & 
Venkatesan 2003). In both analyses, this corresponds to 
the intrahalo timescale for the reincorporation of metals 
created by the first SNe into cold starforming gas clumps 
up to metallicities roughly corresponding to Z cr in primor- 
dial galaxies. The upper limit of 100 Myr corresponds to 
the inter-halo enrichment timescale over which metals are 
transported to and pollute neighboring halos at z ~ 10-20 
(Tumlinson et al. 2004; Bromm et al. 2003). 

Such estimates are consistent with the 1-50 Myr 
timescales needed for the instability in SN shells to set 
in for second-generation star formation (Salvaterra et al. 
2004). Stated another way, this timescale is a requirement 
for the model in this work, where the dust generated in the 
very first SNe, presumably on timescales of order 2-5 Myr, 
is exposed to the hot radiation field from a zero-metallicity 
stellar cluster which is assumed to be within the SNR. We 
focus on the effects of an individual SN in a stellar cluster 
and assume for simplicity that it occurs at the center of its 
host galaxy. A Pop III cluster is likely to form with only 
a few stars in the densest part of the galaxy, and the mas- 
sive stars in the Pop III IMF likely have not migrated too 
far from their birthsites at z 10. The specific location 
may not matter much for this work as we are interested in 
relatively local rather than intergalactic transport. 

A related issue is the spatial overlap and relative centers 
of the SNR and the background radiation field. In the sce- 
nario we propose here, the net radiation pressure within 
the SNR is important, and we assume that the radiation 

4 Selecting merely on the basis of C-cnhanccmcnt is unlikely to narrow down the parent process, as many C-enhanced EMP stars show selective 
enhancements of r- and s-process, as well as other elements. As many authors have discussed, this great diversity of abundance ratios may 
indicate an equal diversity of contributing mechanisms. 
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field of the cluster is not dominated by sources that are 
spherically symmetric external to the SNR (so as to can- 
cel the field in the SNR interior). The Jeans length in 
primordial gas is of order 1 pc, a distance that the radii 
of SN shells in our calculations exceed early in the SNR 
evolution (as shown below). Therefore, we assume approx- 
imate spherical symmetry, and that the Pop III cluster is 
not strongly off-center relative to the SNR. 

The first SNe are required in this model to occur in ha- 
los whose virial temperature is at least 10 3 ~ 4 K, for two 
reasons. First, this is a criterion in the model of triggered 
star formation that we use here (Salvaterra et al. 2004). 
Second, the SN kinetic energies considered here (10 51 -10 53 
erg) do not exceed the binding energies of such galaxies, 
ensuring that we are dealing with subgalactic phenom- 
ena rather than dust transport to the IGM. A single 10 53 
erg SN can expel gas efficiently from low-mass miniha- 
los of virial temperatures 100 K (corresponding masses of 
~ 1O 6 M0); even such small halos can, however, partly 
survive a 10 51 erg SN (Bromm et al. 2003). 

We assume that the dust from the first SNe is created 
on timescales of hundreds of days after the SN event (To- 
dini & Ferrara 2001; Nozawa et al. 2003; Schneider et al. 
2004) and solve for the density and temperature of the 
SNR as detailed below. In young remnants, the dust grains 
presumably form in cooling, metal-rich ejecta, which are 
slowly decelerated by interactions with gas in the SNR in- 
terior and by the reverse shocks. After the dust grains 
form through nucleation of ejecta material, they are sub- 
ject to destruction through thermal sputtering by plasma 
(H, He ions) and by passage of the reverse shock through 
the dense ejecta. We compute the sputtering within the 
SNR (Sedov- Taylor) interior, but do not compute the ef- 
fects of reverse shocks. Some dust grains may survive these 
shocks, if a sufficient number reside in dense cold clumps. 
This may be indicated by observations of the Cas A SNR 
(Greidanus & Strom 1991; Hines et al. 2004), although 
whether the detected dust belongs to the SNR or to cir- 
cumstellar material is controversial (Krause et al. 2004). 

Those dust grains that survive are radiatively driven 
by the Pop III radiation field, and transported to varying 
distances as detailed below. We do not model the effects 
of the inhomogeneous density and radiation field arising 
from multiple stellar clusters, or the nonuniform density 
and velocity structure within the SNR. We assume that 
the initial velocities of all the dust grains are of order 100 
km s _1 , consistent with observations of the Cas A SNR 
(Fesen et al. 1987). We discuss the effects of varying the 
initial velocity of the dust grains in more detail in the next 
section. 

For metal-free Type II SNe (Todini & Ferrara 2001), 
magnetites and ACG (graphites) in general dominate the 
dust mass in the stellar progenitor range 12-35 M , with 
an increasing magnetite yield with rising SN explosion en- 
ergies for 22-35 M©. For the specific case of 22 M©, the 
typical grain sizes are about 0.001 fim for magnetites, 
slightly lower values for silicates, and about 300 A for 
graphites. In contrast, the dust mass from PISNe is 
strongly dominated by silicates, with a small nearly con- 
stant yield of ACG and a rapidly rising yield of magnetites 
with increasing stellar mass over 140-260 M©. In addition, 
a 149 (250) M© parent star generates post-SN dust grains 



of characteristic sizes of 0.001 (0.01) /J,m, 0.01-0.1 (0.001- 
0.01) /<m, and 0.01-0.1 (0.1) ^m for magnetites, silicates 
and graphites. These trends of grain compound and size 
with stellar mass from Nozawa et al. (2003) and Schnei- 
der et al. (2004) use predictions of element nucleosynthe- 
sis from current hydrodynamical SN models. These re- 
sults combined with calculations of the chemistry, tem- 
perature and density evolution of the SN ejecta determine 
which and when dust compounds can condense after the 
SN explosion, and their subsequent growth. The temper- 
ature affects the order of grain formation which in turn 
influences their characteristic sizes. In the PISNe mass 
range, graphites condense first (at higher ejecta tempera- 
tures), and therefore tend to be larger in size relative to 
silicates and magnetites which condense later. Thus, not 
all dust grains are created equal: there is a strong de- 
pendence of grain size and net yield on the star's mass, a 
point whose importance to the EMP stellar abundances is 
demonstrated below. 

3.2. Transport of Dust Grains 

The equation of motion of nonrotating spherical grains, 
of mass md = 4^p s a 3 (with density p g and radius a), un- 
der the effect of radiation pressure, gravity and gas drag 
is given by, 

dv 



m d-7- = F r 
at 



F„-F, 



drag ■ 



(i) 



Here the force F r due to radiation pressure for a given 
source of luminosity L, acting on a grain at a distance r, 
is, 

L 



F r = 



4nr 2 



~7TQ Qpr, • 



(2) 



where Q pr = Q abs + Q sca [l ~ (cos#)]. The averaging is 
done over the spectral energy distribution of the radiation 
field. We assume a primordial composition for the gas. 

Qabs and Q sca are the absorption and scattering coef- 
ficient and (cos 9) characterizes the scattering property 
of the grain. The values of Qabs, Qsca, and {cos 9) for 
graphites and silicates are tabulated in Draine (1985). The 
values of Q SC a are typically much smaller than Q a bs and for 
our calculation wc have used Q sca — 0, which makes the 
estimate of the radiation pressure on grains conservative. 
For the values of Q a bs for silicates, we use the fit provided 
by Ferrara & Dettmar (1994), and for graphites, we use 
the fits provided by Nath et al. (1999), for different ranges 
of incident photon energy. The radiation field is assumed 
to be Planckian with a radiation temperature T* . 

We have performed similar calculations for magnetites 
(Fea04). This compound could be produced in primordial 
SNe (Nozawa et al. 2003; Schneider et al. 2004), subject 
to the uncertainties in the explosion mechanism and mass 
cut of SNe associated with the first stars. For this work, 
the role of magnetites is particularly important to con- 
sider, as the differential transport of graphite and silicate 
dust grains relative to those in iron compounds must be 
accomplished if we wish to explain the origin of the ultra- 
iron-poor EMP stars that have highly enhanced C, O, Mg 
and Si relative to solar values. For magnetites, we have 
used the appropriate values of the absorption coefficients 
calculated from Mie scattering theory (S. Bianchi 2005, 
private communication) . 
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The gravitational force on the grain is calculated as- 
suming a NFW profile (Navarro et al. 1997), assuming the 
source of radiation to be at the cluster center, as, 



charge) , 



GM tot {r)m d 



(3) 



where M tot (r) is the total mass inside the radius r. Fol- 
lowing Komatsu & Seljak (2001), one can define a charac- 
teristic radius r s — r v i r /c for an object of total mass M 
collapsing at a redshift z, and for a concentration param- 
eter c, where, 



M 



(4V3)A c (z)fi m/9c (z) 



(4) 



and A c (z) is the overdensity and p c (z) = 3H 2 /8irG is the 
critical density of the universe. The overdensity A c (z) « 
(18tt 2 + 82x - 39x 2 )/n(z), with x = Cl(z) - 1 and Cl(z) 
as the ratio of mean matter density to critical density at 
redshift z (Bullock ct al. 2001). The density at this char- 
acteristic radius r s is given by, p s — [c 3 M / (47rr 3 ir ,m(c))] 
where m(c) = ln(l + c) — r/(l + c). The total mass within 
a radius r can then be written as, 



M tot (r) =4TT Ps r 3 s m(r/r s ). 



(5) 



To compute the drag force, we first calculate the charge 
on the grain as detailed further in this section, and use 
equation (4) from Drainc & Salpeter (1979a). We refer 
the reader to this paper for greater detail on this calcu- 
lation. The drag force on the dust grain has two com- 
ponents: collisional drag caused by the physical collision 
of grains with H and He atoms, and plasma drag arising 
from the long-range Coulomb forces associated with the 
grain charge. Additional effects for the grain charge that 
arise from the grain's motion (Shull 1978) are relevant only 
when the grain's velocity exceeds the thermal velocity of 
protons, and are not considered here. 

Grains are charged by the photoelectric effect in the 
presence of the radiation field, and collisions with electrons 
and protons. Coulomb drag is the dominant process in this 
calculation and strongly effects the motion and evolution 
of the grains. We do not include the effects of magnetic 
fields in primordial galaxies in this calculation, in the ab- 
sence of a compelling theory as to such a field's structure 
and magnitude and the role of galactic amplification pro- 
cesses. The charging from the photoelectric effect depends 
on Qabs , the absorption coefficient of the grain (see above) 
and the radiation field J„. For a thermal gas with elec- 
tron and proton densities n e ,n p and temperature T, the 
photoelectric current is given by, 



Jph 



f 



Qab s (u,v)yJ^-j- dv , (6) 



where v m i n = (w + eU)/h corresponds to the minimum 
energy for which electrons can escape from the surface, w 
is the workfunction, and y v is the normalized photo-yield. 
We use the fits for y v provided by Bakes & Tielens (1994). 

Following Draine & Salpeter (1979b), we write the cur- 
rents due to collisions with protons, as (with x = eU/kT, 
U = Z gI e/a being the grain potential, and Z gT the grain 



J 



kT \i/2 



and with electrons, as, 



Je 



kT \i/2 



(7) 



(8) 



\2%m e 

These currents are enhanced by Coulomb interactions by 
a factor g(x). For electrons, g e {x) = cxp(a;), for x < 0, 
and g e {x) = 1 + x for x > 0, and for protons, g p (x) = 
exp(— Z gI x) for x > 0, and g p (x) = 1 — Z gI x for x < 
(Spitzcr 1978). 

The charging time scales are much smaller than the 
transport time scale of the grains in our calculations, and 
therefore the grain charge can be calculated from the equi- 
librium condition at each radii : J p h + J p = J e . 

In addition to the equation of motion (equation 1), we 
also calculate the evolution in grain sizes due to sputter- 
ing. Grains are considered to be destroyed when they are 
reduced to 10 -3 of their original size. For the case of mag- 
netite grains, there are no published sputtering yields in 
the literature, and a detailed calculation of this quantity is 
beyond the scope of this work. Bianchi & Ferrara (2005) 
derived that graphite and olivine [(MgjFe^SiCU] grains 
have similar sputtering rates as long as the grain veloci- 
ties are at least 100 km s _1 , as we nhave assumed here. 
The theoretical expectation would be that the values for 
magnetite sputtering lie between those for silicates, which 
have similar binding energies per atom, and for metal- 
lic Fe (Draine 2005, private communication). The latter 
may be preferred if oxygen is preferentially sputtered from 
magnetite grains. For these reasons, we approximate the 
magnetite sputtering rate as that for pure Fe from Tielens 
et al. (1994); for consistency, we use the sputtering rates 
from these authors' calcuations for graphites and silicates 
as well (see also equation 3 of Scab & Shull 1983). 

4. RESULTS 

We present the main findings of our paper in this section 
in four parts in the following order: (1) the role played by 
a grain's area to mass ratio in determining its response to 
a radiation field; (2) solving for the evolution of the densi- 
ties and temperatures within primordial SNRs which will 
provide the background conditions for grain transport; (3) 
using the results from (2) and earlier sections to derive the 
efficient transport of graphite and silicate grains relative to 
magnetites within SNRs, with the associated implications 
for the abundance ratios of Fe-poor EMP stars and the 
first-stars IMF; and (4) quantifying the impact of varying 
a few of our model assumptions on the results in this work. 
At the end of each of the first three subsections, we present 
a brief summary of its findings in the broader context of 
the paper. 

4.1. Radiation Pressure on Dust Grains 

The dynamics of dust grains are determined by a com- 
plex interplay between gas density, temperature, and grain 
properties like size, mass, density, and charging rate. The 
relevant considerations for this problem are: (1) the ratio 
of the grain's effective area to mass, which is important for 
grain dynamics under radiation pressure and gravity; and 



6 



Venkatesan, Nath, & Shull 



(2) the relative values of the outward force due to radia- 
tion and gas drag, another important input to the grain's 
motion. 

The effective area, also frequently referred to as the 
absorption cross section, is determined by the geometrical 
area (na 2 ) and the absorption coefficient Q a bs- We plot 
the ratio of the effective area to the mass of the grains 
(in units of cm 2 g _1 ) as functions of the wavelength of 
the incident radiation in Figure 2 for graphites, silicates, 
and magnetite grains. For each species, we show two cases 
with grain radii a = 0.01 /im and a = 0.01 /im. We should 
note that for a Planckian radiation field with effective 
temperature T e ff = 5 x 10 4 K, the peak of the spectrum is 
at A ~ he/ (2.82 kT e g) ~ 0.1 /im. Such a hot radiation field 
is easily achieved in the atmospheres of metal-free stars of 
masses exceeding about ten solar masses (Tumlinson et al. 
2003; Bromm et al. 2001b). The curves show that smaller 
grains in general have larger ratios of area to mass, and 
graphites have the largest ratio, followed by silicates and 
then magnetites, a direct consequence of their respective 
grain densities of 2.2, 3.2, and 5.2 g cm -3 . 




-10 12 
Log X (fJ-tn) 

Fig. 2. — The ratio of effective area to grain mass is plotted as a 
function of wavelength of the incident radiation, for graphites (solid 
lines), silicates (dotted-dashed), and magnetite (dashed) grains. In 
each case, thick and thin curves denote grains of size a = 0.001 fim 
and a = 0.01 fim. 

We therefore anticipate, in the context of SNRs result- 
ing from Pop III SNe in primordial galaxies, that graphite 
grains will be affected by radiation pressure more than sili- 
cate and magnetite grains, and small grains will experience 
greater radiation pressure than large grains. 

4.2. SNR Evolution 

In order to calculate the transport of dust grains with 
the equations described in the previous section, we need 
to assume approximate background gas temperatures and 
number densities within the SNR. Although we have not 
modeled the detailed density and velocity structure in this 
region, we assume parameter values that encapsulate the 
range of physical conditions expected within the SNR in 
the SN-induced star formation scenario. We show this 
in Figure 3, where the evolution of the density and tem- 
perature with interior radius in the SNR (and not the 
SN shell radius) is plotted for times up to 10 Myr. We 



use Sedov- Taylor self-similar solutions and the analytic 
approximations of Kahn (1975) and Petruk (2000) for the 
temperature and density profiles in Figure 3, assuming a 
typical HN/PISN explosion energy £?sn = 10 52 erg and an 
ambient gas density of 10 cm~ 3 . These parameters repre- 
sent the typical conditions under which the fragmentation 
instability can occur in SN shells in primordial galaxies 
(Salvaterra et al. 2004), although densities of 10 cm~ 3 
may only be achieved in the cold star-forming cores of 
such galaxies. 




Log rsNR (pc) 

Fig. 3. — The evolution of the particle density (solid lines; left 
y-axis) and temperature (dashed-dotted lines; right j/-axis) in the 
SNR interior with radius in pc, for an ambient density of 10 cm -3 
and -Esn = 10 52 erg, and for an assumed Sedov- Taylor self-similar 
solution. Within each set of lines, the curves from left to right are 
for t = 10 4 , 10 s , 10 6 and 10 7 yr; these would move to the right for 
lower ambient densities (see text). 

The timescales on which SN shells form can be estimated 
from the equations describing the Sedov- Taylor phase of 
SNRs (Shull & Silk 1979). For the above parameters, the 
shell radius, velocity, and post-shock temperature in the 
shell scale with E52 (SN energies of 10 52 erg), ambient 
density n\o (in units of 10 cm -3 ) and time £4 (in units of 
10 4 yr) as: 

R s = [12.7 V c) E»i n^-h^ (9) 
V s = (439 km s- 1 ) EHn^% °- 6 (10) 
T s = (2.58 x 10 6 K) E^n^-% 1 - 2 (11) 

Applying the cooling criterion for the onset of radiative 
shell-formation with free-free cooling from H + and He +2 
ions leads to shell formation at times of: 

Uell ~ (2.1 x 10 4 yr) E^n^ 4 (12) 

Thus, the SN shell has begun to form at about 10 4 
yr, and the temperature has dropped sufficiently to lower 
sputtering rates within the shell. Figure 3 shows that, un- 
der these conditions and for the SN energies considered 
here, the typical size of SN shells is on the order of tens to 
at most hundreds of pc over timescales of 10 Myr. We also 
investigated cases with ambient densities of 1 cm -3 , where 
the lines in Figure 3 simply scale towards the right. At a 
given time and SN shell radius, the density within the SNR 
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is almost two orders of magnitude lower in comparison to 
the 10 cm -3 case. 

Figure 3 also reveals that at early times (< 10 5 yr) when 
the SN shell has advanced to about a few to 10 pc, typical 
densities and temperatures in the SNR interior are respec- 
tively on the order of 0.1 cm~ 3 and a few times 10 8 K. 
Since the dust that survives post-SN processes such as the 
reverse shock is likely to be found in cooler gas, we conser- 
vatively assume ambient temperatures of 3 x 10 7 K (the 
value at about 6 pc at 10 4 yr), and an ambient density of 
0.1 cm' 3 within the SNR. 

In summary, for a typical Pop III SN with E§n = 10 52 
erg, a SN shell begins to form in primordial galaxies at 
~ 10 4 yr after the initial explosion, and reaches distances 
of ~ 10-100 pc in 10 Myr. The SNR interior remains hot 
on these timescales, and we assume background gas con- 
ditions of 3 x 10 7 K and 0.1 cm -3 in solving for dust grain 
transport within the SNR. 

4.3. Grain Dynamics 

We use the results of Schneider et al. (2004) and To- 
dini & Ferrara (2001) for the formation of dust grains in 
SNe in primeval galaxies in order to assume the sizes of 
grains of different composition in our calculations. They 
found, as discussed earlier in this work, that for SNe 
with a progenitor mass of 149 Mq, the size distribution 
of magnetite grains peaked at ~ 0.001 fim, and those of 
graphites and silicates at ~ 0.04 /im. For stellar masses 
of 22 Mq and 250 Mq, the characteristic grain sizes for 
magnetites/silicates/graphites are respectively of order 
0.001/0.001/0.03 jum and 0.01/0.01/0.2 fj,m. These three 
masses span the range of SN progenitors considered here. 
We refer the reader to these papers and to Nozawa et al. 
(2003) for the detailed grain size distribution for various 
stellar masses. 

1.5 E 1 ' ' 1 ' 1 



1.0 




-0.5 E , , , , , , , 3 

-2.4 -2.2 -2.0 -1.8 -1.6 -1.4 -1.2 -1.0 
Log t(Myr) 

Fig. 4. — The positions of dust grains within SNRs in primor- 
dial galaxies are shown as a function of time as they are radiativcly 
transported outward from the center of the SNR. Top, middle and 
bottom panels represent stellar masses of 22 Mq, 149 Mq, and 
250 Mq, spanning the range of SN progenitors considered in this 
work. Solid, dotted and dashed curves refer to graphites, silicates 
and magnetite grains. The initial grain sizes for the compounds in 
each panel are detailed in the text. The two dashed-dottcd lines in 
each panel describe the adiabatic evolution of the SN shell radius, 
for -EgN = 10 52 erg and and ambient gas number densities of 1 cm -3 
(upper line) and 10 cm -3 (lower line). See text for discussion. 



We show in Figure 4 the position of dust grains under 
the combined forces of radiation, gas drag, gravity, and 
sputtering, using the above values for the initial sizes of 
grains of the relevant composition for the three represen- 
tative stellar mass cases. We recall our model assumptions 
discussed earlier: that all grains are released at r = with 
an initial speed of 100 km s _1 , that the luminosity of the 
central source is L — 10 6 Lq with a Planckian spectrum of 
effective temperature T c g = 5 x 10 4 K, and that the tem- 
perature and gas density within the SNR are uniform at 
3 x 10 7 K and 10 _1 cm -3 respectively. We do not use the 
temperatures and densities from the self-similar solutions 
shown in Figure 3 for two reasons. First, strong departures 
from a roughly constant product of density and temper- 
ature within the SNR occur only close to the SN shell. 
Second, the self-similar profile would indicate extremely 
high temperatures (> 10 10 K) near the SNR center, which 
would destroy all dust grains promptly. Such a profile is an 
idealization in the very central regions, and it may not be 
a realistic representation of the physical conditions there. 

Figure 4 shows the radiative transport of graphites, sil- 
icates and magnetites, as well as two additional curves 
that represent the adiabatic evolution of the SN shell for 
an explosion energy of 10 52 erg and ambient densities of 1 
cm -3 (upper curve) and 10 cm -3 (lower curve). Clearly, 
the SN shell advances to smaller distances for higher am- 
bient densities, which could lead to more effective pile-up 
of grains. The SN shell positions are computed using the 
Sedov- Taylor formalism in the appendices of Shull & Silk 
(1979) for the first 10 2 yr, and thereafter by solving the 
differential equations (2) and (3) in Salvaterra et al. (2004) 
without the radiative cooling term in eqn. 3. We note that 
the latter treatment accounts for the mass swept up by the 
shell, which is substantial by 10 4 yr and considerably slows 
the forward motion of the shell. Hence the shell radii val- 
ues in Figure 4 are significantly smaller than those from 
a simple Sedov- Taylor treatment for times exceeding 10 4 
yr, e.g., 11 versus 19 pc at 10 4 yr for ambient densities of 
1 cm~ 3 . 

The curves in Figure 4 show that dust grains of sufficient 
size are driven by radiation equally effectively irrespective 
of grain type, and that graphites, silicates, and magnetites 
may pile up in the SN shell well before 1 Myr. Although 
smaller grains experience radiation pressure more strongly 
(Figure 2), they are destroyed on faster timescales owing 
to greater relative loss of size. This is seen clearly in Figure 
4, where graphites consistently reach the SN shell owing 
to their relatively large size for all the displayed stellar 
masses, whereas the small magnetite grains are destroyed 
at the high SNR interior temperatures before arriving at 
the shell. As the temperature increases, the sputtering 
rates rise and compete with the decreasing Coulomb drag. 
Sputtering affects graphite and silicate grains as well, but 
given their typical grain sizes in the figure, significantly 
higher densities and/or temperatures (> 7 x 10 7 K) are 
required to boost sputtering and prevent shell pile-up of 
these grains. This can also be achieved by having the SNe 
explode into an ISM of lower ambient densities than we 
have assumed here, in which case the SN shell would ad- 
vance to larger distances, leading to less effective pile-up 
of dust grains. Transport through SNR interior temper- 
atures that exceed ~ few xlO 7 K will destroy grains of 
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similar sizes at roughly the same radii, and will not cause 
element segregation. 

In the context of EMP stellar abundances, where silicate 
elements and carbon were often enhanced with respect to 
iron, the 149 Mq case appears at first glance to offer the 
best explanation for reasonable parameter choices. The 22 
Mq case favors graphites alone, whereas a SN from a 250 
Mq star creates larger grains in general, leading to the 
eventual transport of graphites, silicates and magnetites 
to the SN shell. The magnetite grains from 22 Mq and 
149 Mq SNe are considerably smaller (~ 0.001 /mi) and 
are sputtered away before they can arrive at the SN shell. 

To fully reconcile these results with the current EMP 
data presented in §2, there is an additional important fac- 
tor to consider: the relative masses in these dust grains 
formed in primordial SNe. The survival and deposition 
of any compound must be considered along with its net 
yield. As discussed earlier, magnetites and graphites dom- 
inate the dust yield from 12-35 Mq SNe, while PISNe from 
the stellar mass range 140-260 Mq create mostly silicates, 
typically 1-2 orders of magnitude more by final mass than 
graphites or magnetites. When these trends are consid- 
ered with those in Figure 4, we conclude that the mass 
range ~ 10-150 Mq best explains the enhancement of car- 
bon and silicate elements in Fe-poor EMP stars, with the 
lower end of this range being most appropriate for C-rich 
Fe-poor EMP stars. High-mass PISNe, as seen in the 250 
Mq case, would predict a strong enhancement of silicates 
and iron relative to carbon, which is not matched by the 
data. 

These results remain approximately true when com- 
pared with the predictions of dust grain sizes and mass 
yields from other works such as Nozawa et al. (2003). 
These authors predict larger grain sizes and yields for 
silicates from a 25 Mq star than we have assumed here 
for a 22 Mq star from Todini & Ferrara (2001), and a 
mass output of graphites that depends strongly on the de- 
gree of mixing in the SN ejecta. The trends from Nozawa 
et al. (2003) would imply that more silicates from metal- 
free Type II SNe would reach the SN shell, and that the 
transport of graphites from 10-260 Mq stars would be 
more variable, depending on the physical conditions in in- 
dividual SNe. These outcomes will only strengthen the 
connection to the EMP stellar data, where the elements 
Mg, Si and O are typically overabundant relative to Fe at 
low [Fe/H], while this is true for C in only a fraction of 
the same EMP stars. 

To summarize, current models indicate that primordial 
SNe create graphite and silicate grains of larger size and, 
for some stellar masses, in greater quantities than mag- 
netites. This directly results in the more efficient transport 
of graphites and silicates through the SNR interior to the 
SN shells, where they could provide cooling that triggers 
second-generation star formation. A comparison with the 
abundance trends in Fe-poor EMP stars indicates a first- 
stars IMF in the mass range 10-150 Mq, if dust transport 
dominates the element ratios in these stars. 

4.4. Model Assumptions and Variations 

We now address a few remaining issues related to our 
model assumptions, starting with our estimate of the mag- 
netite sputtering rate. The calculation of this quantity has 



two input parameters - the density and sputtering rate of 
a specific grain species. We recall that currently there ex- 
ists no sputtering data for magnetites; the rates for either 
silicates (which has a similar binding energy per atom) or 
for pure metallic Fe can be used as an approximation. The 
magnetite grain density may be set to its measured value 
of 5.2 g cm~ 3 , or to that of pure Fe (7.9 g cm -3 ). There 
are therefore four possible combinations of these two pa- 
rameters in calculating the sputtering of magnetite (recall 
that for magnetites in the above figures, we have assumed 
grain densities of 5.2 g cm" 3 and a pure Fe sputtering 
rate). 
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Fig. 5. — The radiative transport of magnetite grains is shown as 
a function of time with varying parameter values for grain density 
and sputtering rates. Dotted and dashed lines represent sputtering 
rates corresponding to those for silicates and pure Fe respectively; 
within these cases, upper and lower curves correspond to densities of 
7.9 g cm -3 and 5.2 g cm -3 (see text). The case assumed through- 
out the paper is marked with an arrow and is intermediate to the 
other cases. The initial size of all magnetite grains is assumed to be 
0.001 nm. 

We show in Figure 5 the time-dependent position of 
magnetite grains for these four parameter combinations 
for grain sizes of 0.001 /im. The curves which use silicate 
sputtering rates in general lie below and terminate ear- 
lier than those with pure Fe sputtering rates, as expected 
from the higher values of the first quantity. A perhaps less 
obvious result in Figure 5 is the larger distance to which 
the grains with higher density are transported. This re- 
veals the complex nature of the problem of dust transport; 
as we emphasized earlier, individual trends of grain mass, 
size, and density do not combine predictably in a given 
problem. In this case, the difference between the trajec- 
tories of light (less dense) and heavy (more dense) parti- 
cles is not simply extrapolated, as they are subjected to 
a position-dependent radiation field. One may have ex- 
pected that heavier particles always have lower velocities 
than lighter particles. We calculated velocities for each of 
the cases in Figure 5 and find that the lighter particles do 
indeed move faster initially, but slow down more rapidly 
than do the grains with higher densities. We attribute this 
behavior to the effect of the decreasing values of the driv- 
ing radiation field with distance. Another factor may be 
the increase of grain charge with decreasing grain density. 
This trend appears to hold for the physical conditions in 
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the problem explored in this paper, where the charge in- 
creases from magnetites to silicates to graphites, leading 
to the greater slowing of lighter dust grains from gas drag. 
This is consistent with the results of Bianchi & Ferrara 
(2005) in the context of dust transport to the IGM. We 
emphasize that this anti-correlation between charge and 
grain density may be true only in the specific environment 
considered here, and may not apply in other astrophysical 
situations. Our main point which we wish to demonstrate 
through Figure 5 is that the curve corresponding to grain 
density 5.2 g cm -3 and a pure Fe sputtering rate lies inter- 
mediate to the other cases in the figure, and justifies our 
earlier assumption of these values for magnetite grains. 

Last, we have explored variations of the initial velocities 
of the dust grains. Lowering this quantity from our as- 
sumed value of 100 km s _1 to 10 km s _1 does not change 
our results appreciably, because the speed to which the 
grains are accelerated by radiation pressure is itself of or- 
der 100 km s _1 . However, higher initial velocities of order 
1000 km s _1 generally result in the grains being trans- 
ported to greater radii. Depending on the role played by 
grain charge and gas drag, the grains may driven out of 
the galaxy (see earlier references on dust transport to the 
IGM) , and become less relevant for the problem of seeding 
metals for second-generation star formation in early halos. 

5. CONCLUSIONS 

We have investigated the radiative transport of dust 
within SNRs in primordial galaxies. We find that the dif- 
ferential transport of the primary dust compounds from 
the first SNe could be important for the problem of which 
metals seed the second generation of stars that form in the 
cooling shells of such SNe. The survivors of such second- 
generation star formation may be detected as the EMP 
stars in the Galactic halo, whose metal abundance ratios 
provide clues to the chemical environment in which they 
formed. Calculations of the motion of dust grains within 
SNRs under the combined forces of the hot first-stars radi- 
ation field, gas drag, gravity and sputtering reveal a com- 
plex interplay between grain size, charge, mass, density, 
and species. Our model requires that Pop III star forma- 
tion can continue for at least a few to 10 Myr in early 
halos, the lower limit on the sum of the lifetimes of mas- 
sive stars and of the timescale for dust grains to cross the 
SN shells (typically well before 10 5 yr). In the context of 
this problem, our main findings are: 

1. The role played by radiation pressure in the grain 
dynamics of different species is strongly influenced 
by the ratio of the grains' area to mass, and de- 
creases from graphites to silicates to magnetites. 
Within a given species, small grains are driven more 
effectively by radiation than are large grains. 

2. We find that grain charge increases from magnetites 
to silicates to graphites, leading to greater Coulomb 
drag and loss of velocity for graphites and silicates 
relative to magnetites. However, a competing effect 
is the faster destruction of small grains, despite the 
stronger radiation pressure experienced by grains 
of decreasing size. The net result is that the gener- 
ally larger graphite and silicate grains from the first 



SNe are more effectively driven and accumulate in 
the SN shells more than small magnetite grains. 

3. If a dust-transport scenario in primordial SNRs 
were to account for the segregation of elements ob- 
served in C- and silicate-rich Fe-poor EMP stars in 
the Galactic halo, a metal-free stellar IMF spanning 
~ 10-150 Mq offers the best explanation. 

4. Higher temperatures and densities than those we 
have considered here will lead to increased sputter- 
ing and destruction of most grains, and a loss of se- 
lective grain transport. At lower ambient ISM den- 
sities, the SN shell would advance further, reducing 
the possibility of grains of any species reaching the 
SN shell and mixing with star-forming gas. 

A detailed prediction of the contribution of dust to the 
metallicity of second-generation stars is beyond the scope 
of this paper, and is likely best studied through numer- 
ical simulations of primordial star-forming sites such as 
those cited earlier in this work. Clearly, dust transport 
can enrich such sites in elements that constitute the grains 
reaching the SN shell. A 200 M Q star can produce up to 
60 Mq of dust (Schneider et al. 2004); if this is trans- 
ported completely to the SN shell at a radius of 10 pc 
with an ambient density of 10 cm~ 3 and diluted into the 
swept up mass of ~ 10 3 M Q , then supersolar mctallicities 
of about 3 Zq in the shell can result. This is only an es- 
timate of the average value. The real metallicity may be 
much higher or lower in cold star-forming clumps in the 
SN shell, depending on the efficiency of mixing processes 
and the role played by some of the transported elements in 
providing cooling pathways for primordial gas. Although 
it is beyond the scope of this work to fully compute how 
and when the transported dust grains get mixed into the 
SN shell gas, we may estimate it from the gravitational 
instability or free-fall timescale. This quantity is approxi- 
mately {Gp)^ 1 / 2 , and equals 10 Myr for hydrogen number 
densities of ~ 10 2 cm~ 3 . Therefore, if there is to be no 
significant delay from dust grain mixing prior to the on- 
set of second-generation star formation, the gas in the SN 
shell has to exceed number densities of about a few x 10 3 
cm~ 3 . 

There remain sizeable uncertainties in the modeling 
of dust transport in astrophysical environments, such as 
those detailed above, and factors related to the creation 
of dust in primordial SNe such as the role of reverse 
shocks. Clearly, some dust does survive in these envi- 
ronments, as inferred from the recent detection of dust 
inaz~ 6.2 QSO whose extinction curve indicates a SN 
origin (Maiolino et al. 2004). An additional point that 
we have not accounted for is the binarity fraction of the 
first stars, which can affect the mass loss from companion 
stars, and reverse shock phenomena related to circumstel- 
lar material. Future spectroscopic analyses of Sloan-DSS 
and Hamburg/ESO survey data on EMP stars and ob- 
servations of a variety of Galactic SNRs will help to ad- 
dress these important issues and ultimately constrain the 
formation sites and conditions of cosmological first- and 
second-generation star formation. 
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